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Fluorescence data can be rapidly acquired in the form of an emission-excitation matrix (EEM) using a novel 
fluorometer called a video fluorometer (VF). An EEM array of 4096 data points composed of fluorescence intensity 
measured at 64 different emission wavelengths and excited at 64 different excitation wavelengths can be acquired in less 
than one second. The time-limiting factor in using this information for analytical measurement is the interpretation step. 
Consequently, sophisticated computer algorithms must be developed to aid in interpretation of such large data sets. For 
"r" number of components, the EEM data matrix, M, can be conveniently represented as 

r 

i = l 

where x(0 and y(0 r are the observed excitation and emission spectra of the i^ component and a t is a concentration 
dependent parameter. Such a data matrix is readily interpreted using linear algebraic procedures. 

Recently a new instrument has been described which rapidly acquires fluorescence detected circular dichroism 
(FDCD) data for chiral fluorophores as a function of multiple excitation and emission wavelengths. The FDCD matrix 
is similar in form to EEM data. However, since the FDCD matrix may have legitimate negative entries while the EEM 
is theoretically non-negative, different assumptions are required. This paper will describe the mathematical algorithms 
developed in this laboratory for the interpretation of the EEM in various forms. Particular emphasis will be placed on 
linear algebraic and two-dimensional Fourier Transform procedures. 
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Introduction 

Advances in computer technology and developments in 
multiparametric detection devices have had a profound ef- 
fect on developments in chemical analysis.'' These develop- 
ments have made it possible to expand the applicability of 
several analytical methods to more complex systems. Multi- 
component analysis by fluorescence spectroscopy is one 
example of such a method. 
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In the conventional fluorescence experiment, the sample 
solution is irradiated with monochromatic light which pro- 
duces molecules in the excited state when absorbed. As 
these molecules return to the ground state, light with a 
characteristic wavelength distribution is emitted. This distri- 
bution of the emitted light is known as the fluorescence 
emission spectrum. The fluorescence excitation spectrum is 
the fluorescence intensity as a function of absorption wave- 
length. At low absorbance (<0.01), the intensity of the 
fluorescence emitted at a given wavelength is directly pro- 
portional to the amount of light absorbed and therefore to the 
concentration of the analyte in the sample solution. 

These characteristic spectra make fluorescence spec- 
troscopy inherently more selective than absorption methods 
and provide qualitative as well as quantitative measure- 
ments. For example, when the excitation spectrum of the 
analyte only partially overlaps the excitation spectrum of the 



487 



other components in the sample, the sample can he irradi- 
ated with light of a wavelength which is only absorbed by 
the analylc. Since only the analyte absorbs the incident 
lighi. fluorescence svjj] on)}' be emitted by the analyte. in Ihe 
absence of synergistic effects such as energy transfer. This 
technique is called selective excitation 1 1 1' and can be used 
qualitatively by acquiring the complete emission spectrum 
ol' the analyte or quanliiativcJy by determining the analytc 
concentration using the emission of a sample of known 
concentration as a standard. 

This approach to muliicomponent analysis can be ex- 
panded hy acquiring fluorescence spectra at several excita- 
tion and emission wavelengths. Acquiring multiple spectra 
with conventional i:is[:xii]icma»jtir is a titv.e-cc-nsuming pro- 
cess, even witli a niicnrproccssor-eorilrtilleu instrument. It 
can require more than one hour to acquire the complete 
emission spectrum of a sample at 64 excitation wavelengths 
using conventional instrumentation. Many samples would 
undergo significant photodecomposition over such a period 
of time. The development of the videufluorometer (VF) 
which rapid! v acquires two dimensional fluorescence data 
has resolved this problem [2\. This instrument uses poly- 
chromatic excitation ;tnd a silicon intensified target vidicon 
detector (television camera) to acquire data in matrix format 
without mechanical scanning. The VF can acquire 64 emis- 
sion spectni generated at 64 excitation wavelengths in less 
than 1 second. Therefore, data processing of this emission- 
excitation matrix (fiEMl is now the time-limiting step in the 
analysis and requires the use of a computer for data reduc- 
tion. 

The rows of the HEM of a pure sample are multiples of 
the emissiun spectrum of that compound, while the columns 
yre multiples of the ctcitiition. spectrum. tndicutina Inn: the 
EEM is of bilinear form. Sucii mm rices are ideally suited for 
data reduction techniques such as factor analysis |3] and 
pattern recognition [4|. 

Aiioihcr technique whose applicability is expanded by 
niiillipanimelrie detection is fluorescence detected circular 
dichroism (IDCDl \5] The FDCD spectrum is the differ- 
ence in the intensity of flurnescence produced by excitation 
with left and right circularly polari/ed light recorded as a 
function of the wavelength of the exciting light. Chiral 
molecules preferentially absorb one form of circularly polar- 
ized light. The FDCD spectrum reflects Ihe structure of the 
chiral Jluorophore in a solution. A two-dimensional rapid 
scanning FDCD spectrometer has been developed which 
measures ihe FDCD at several excitation and emission 
wavelengths [6|. When these data are collected in matrix 
format, they arc also of bilinear form. However, the al- 
gorithms designed lo resolve the spectni of components 
from the HUM cannot be applied to this elliptic ity matrix 
since it may con I air. legitimate negative values. Thij. 



manuscript will discuss the strategies developed for the 
qualitative reduction of muJticomponent MEMS as well as 
the alternative techniques developed for multieomponent 
ellipticity matrices. 

Eigenvector Analysis [7] 

When the imsorbanee of u sample is less than 0.0 1, the 
intensity of the fluorescence, /,-. can be approximated by the 
expression 



/, = 2.303 Infydx- 



(la) 



where /,. is the intensity of the incident radiation. <$>.■ -s the 
rl Lioresccnce quantum efficiency fthe traction of absorbed 
photons emitted as fluorescence), e is the molar extinction 
coefficient, h is the thickness of sample cell and c is the 
concentration of the fluoraphoTe in the sample solution. 
Each element, hi,,-, of the emission-excitation matrix. M. 
represents the fluorescence intensity at wavelength \j that 
wits generated by excitation at wavelength X s . Therefore, 
each of these elements can be generally expressed as 



m ti = 2.3034y/ ( A; )e(X,-h<^ M^j Mw 



(2a) 



where -yi X / 1 reflects the dependence of ! t on the monitored 
emission wavelength and 8(X ; -J is a parameter which incor- 
porates instrumental artifacts like sensitivity and signal col- 
lection geometry. Combining these terms based on excita- 
tion and emission wavelength related variables results in a 
simpler expression: 



nifi=ax l y l 



<3r.i 



where a is a scalar that equals 2.303<j> ; Jk" , -r, is the excitation 
term given by 



x^J^kiHtKi) 



and Yj is the emission term which is expressed 



y ; =^A,)S(X,) 



(4a) 



(5a) 



When the .v,- are properly sequenced, the array of.r, is a 
representation of the excitation spectrum and can be denoted 
x in vector notation. Likewise, when the y are properly 
sequenced, the vector y represents the emission spectrum. 
Since (he emission profile is independent of the exciting 
wavelength, and the excitation profile is independent of the 
monitored emission wavelength, the matrix M can clearly 
be expressed as the vector product of x and y, multiplied by 
a scalar concentration term. i.e. 



'Viauivs ii'i br.iiAitw intlK-.ita liieraaw,: pjfcienjcv 



.Vl=ctjcy 



(oa) 
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[n a sample containing "n" fluorescent compounds, the ma- 
trix M is the sum of the EEMs of the individual components 
provided synergistic effects are negligible. Thus, the n com- 
ponent matrix can be expressed as 



M 



2 «* x *y* r 

k = \ 



(7a) 



A more convenient notation for M is 



M=XY 



(8a) 



in which the columns of the matrix X are the excitation 
spectra, \ k , of the n components, and the rows of Y are the 
emission spectra, y k . of the components. The concentration 
term can be considered to be absorbed into either of the 
matrices. 

Qualitative analysis of the matrix M requires a determina- 
tion of the number of independently emitting compounds, 
i.e., the rank of the matrix, and the set of basis vectors x t 
and y k which are the excitation and emission spectra of the 
components. Eigenanalysis plays a significant role in both 
determinations. Therefore, it is appropriate to preface the 
discussion of rank estimation and spectral resolution with a 
brief presentation of pertinent eigenanalysis principles. 

An eigenvector is defined as any vector x which is a 
solution of the equation 



Ax = \x 



(9a) 



in which \ is a scalar called the eigenvalue. The magnitude 
of the eigenvalue is a consequence of the importance of the 
information reflected in the eigenvector to the data in the 
matrix. When the eigenvalue is large, then the factor repre- 
sented by the eigenvector makes a large contribution to the 
data. If it is small, then the contribution of the factor is 
small. Therefore, when a matrix has a rank n which is 
greater than 1 , it has n eigenvectors and n eigenvalues. 

Since the matrix M is bilinear, the covariance matrices, 
M r M and MM 7 " can be used to generate the eigenvectors. 
The eigenvalues of the covariance matrices are the squares 
of the eigenvalues of M, but apart from this small detail, it 
is more expedient to use the covariance matrices to generate 
the eigenvectors since the covariance matrices are always 
square and symmetric. 

As the preceding discussion indicates, rank estimation of 
an ideal matrix (one that is noise-free) is straightforward, 
simply determined by the number of non-zero eigenvalues. 
However, experimental data matrices are not free of noise. 
They can contain systematic or random errors superimposed 
on the signal. For these matrices, there are several methods 
of rank estimation [3,8], When the signal-to-noise ratio is 
high, it is possible to differentiate the eigenvalues associated 



with fluorescence (primary eigenvectors) from those associ- 
ated with noise (secondary eigenvectors) by a direct com- 
parison of their magnitudes. As the signal-to- noise ratio 
decreases, (his approach becomes more difficult Another 
method which will be described in a later section of this 
manuscript is differentiation of eigenvectors based on their 
frequency distributions, realizing that random noise in spec- 
tra! data usually will have higher frequency than the signal. 
The other part of the analysis of M is the resolution of the 
basis vectors x k and y k from the matrix. An infinite number 
of basis vectors exists for a given matrix, and resolving the 
spectra! vectors from the matrix is usually not possible with- 
out a priori knowledge of the components. However, the 
eigenvectors arc an orthonormal basis for M and arc easily 
generated. Since the eigenvectors reproduce the matrix, M 
can be expressed as 



M-UV 



(10a) 



in which the columns of U are the excitation eigenvectors, 
u k , and the rows of V are the emission eigenvectors, v { . The 
eigenvalues have been absorbed into one of the matrices. 
Since the eigenvectors are an orthonormal basis, ihey often 
contain negative elements even though M is theoretically a 
non-negative matrix Emission and excitation spectra are 
also theoretically non-negative. Therefore, the eigenvectors 
can be transformed to possible spectral vectors by trans- 
forming them to a non-negative basis set. 

To perform this transformation, the values of the matrix 
K and its inverse, K ', which transform U and V to non- 
negative matrices, must be found. This is algebraically 
sound, since it is equivalent to multiplying U and V by the 
identity matrix. This transformation is mathematically ex- 
pressed as 



M-UKK 'V 



(Ha) 



This condition also ensures that the transformed vectors are 
also a basis for M. The transformed excitation vectors are 
the columns of the matrix U', and the transformed emission 
vectors are the rows of the matrix V. These matrices are 
given by the equations 



and 



U'=UK>0 



V' = K-'V^0 



(12a) 



(13a) 



The values of the elements of K and K ' can be found from 
the expressions for the elements of U' and V. In the two 
component case, the elements are given by 
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and 



"21 = k n u 2i +k n it u ^Q , 

v' y = U\k\(k 22VlJ -k l2 v v )^0 , 

v' 2j = U\k\(k n v 2j -k 2lVl j)>0 . 



(Ma) 
(15a) 
(16a) 

(17a) 



and 



It can be assumed that k { \=k 22 = 1 without loss of generality. 
The values for k u and k 2] can be found by solving these 
expressions. It's clear from these expressions why this 
method is only applicable to two component matrices since 
the elements of K" 1 become non-linear for more than two 
components. The boundaries for the values of k t2 and k 2l 
given by these expressions are 



v 2; >0 



Vy 



s £i 2 s= max 



~"2, 
Hi, 



(18a) 



min Vy ^2 I s5 max 
v i} u 2i >0 



(19a) 



Vy>0 



The accuracy of transformations performed with values 
meeting these criteria has been shown to be a function of the 
overlap of the spectra of the components. This is illustrated 
in the ambiguity table in figure 1 which summarizes the 
results of the transformation for the 16 possible spectral 
overlap combinations for a binary mixture. It should be 
noted that in 7 out of the 16 possible cases, at least one 
spectrum of each component is correctly resolved. Close 
inspection of the table shows that in cases where the spectra 
of neither component are enveloped in both dimensions, at 
least one of the spectra of each component is resolved unam- 
biguously. Clearly, if both the emission and excitation spec- 
tra of the two compounds have only partial overlap, all four 
spectra will be resolved unambiguously from their mixture 
matrix. 



Emission Overlap 



com p. 1 . 

comp. 2 



a 
n 

> 
O 



O 

Ul 



mln _m. _ ka _ max _H 
v 2| >0 V 2 , u,,>0 u, 



mln 3. - 
Vi,>0 v,, 

All spectra certain. 



max jm 
u 2 i>0 u 2 . 



v 2| >0 v,. 



mln Ji 
v,.>0 ¥,. 



_ max _ 
u,,>0 Ui, 



k;, > max — 
u 2 ,>0 u 2l 



mln 2lL > 
v 2| >0 v 2| 



mm _^L 
"i,>0 v,. 



> max _i 
u N >0 u,i 



, max _J1 
u 2 i>0 u 2l 



1 spectrum of ea comp. certain 
1 spectrum of ea component 
given by extreme k. 



mln Ji _ k 12 
v 2j >0 Vj, 



mln 12. 



Ui,>0 u„ 



> max — 

u 2 ,>0 u 2l 



1 spectrum of ea camp, certain. 
1 spectrum of ea component 
given by extreme k. 



mln Jl > k ,j 
v 2l >0 v 2| 



mln 3. 

Vi,>0 V,, 



^ max _ 
u„>0 u„ 



u a >Q uj. 



I spectrum of ea. comp. certain, 
t spectrum of ea. component 
given by extreme k. 



mln L > k 12 S max 

v 2 ,>0 v 2l u 1( >0 u n 



mln 2L > k,, 

V, : >0 V, ; 



-Ui, 

> max _ 

u 2i >0 u 2l 



All spectra uncertain. 

Both spectra ol ea component 

given by extreme k's. 



mln _li 
v 2) >0 v 2| 



mln _!l 
v,,>0 v. 



k, 2 = max — 
u,,>0 u,, 



k ;1 > max — 
u 2l >0 U;, 



1 spectrum ol ea comp. certain. 
1 spectrum of ea component 
given by intermediate k. 



mm _i > k 12 > max _ 
v 2| >0 v 2i u,,>0 u,i 



v 2 . - u,i 
mln -J. > k 31 > max 1 



v,,>0 v,. 



Uj,>0 u 2! 



mln ^i > k 12 > ma!t "— 
v 2j >0 v 2| Un >0 u„ 



mln 1 

v,,>0 v,, 



- u,, 

max _ 
u 21 >0 u 2 . 



1 spectrum of ea comp. certain. 
1 spectrum of ea component 
given by extreme k. 



min _L > k,j 

v,>0 v 2 



v,,>0 v,, 



> max ™ 

u„>0 u,, 



max 

u 2l >0 u 2 , 



1 spectrum of ea comp. certain. 
1 spectrum of ea component 
given by intermediate k. 



mln _li > K , 3 > max _ 
v 2| >0 v 2l U,i>0 u 1( 



v 2 . 
mm Jl > K 

v,,>0 v,, 



> max _ 
u 2l >0 u s , 



All spectra uncertain. 

Both spectra of ea component 

given by extreme k's. 



All spectra uncertain. 

Both spectra of ea component 

given by extreme k's. 



All spectra uncertain. 

t spectrum of ea comp. given 

by ext. k, other by int k. 



min _i > k, 2 > max 

v 2| >0 v 2| u„>0 u„ 



min _2 > * > max _J1 
Vi,>0 v., u 2i >0 u 2l 

All spectra uncertain. 

1 spectrum of ea comp given 

by ext. k, other by Int k. 



Figure 1— Ambiguity Table of Two-Component Excitation-Emission Matrix. 



min -^ > k, a > max ^1 
v 2| >0 y 3| u „>o u„ 



ml " 21 a k 2 , > max ^ 
v,|>0 v,, u sl >o u 2 . 

All spectra uncertain. 

Both spectra of ea component 

given by extreme k's. 



Vl 

mln Jl > k „ > 
v 2 ,>Ov 2| 



-uj, 
U,i>0 u,. 



mln ^! > k > max "^ 
1 v. — *' ■■ ^n i. 



v„>0 v,. 



u 2 i>0 u 2 . 



All spectra uncertain. 

1 spectrum of ea comp. given 

by ext. k, other by Int k. 



v i| -u 2 , 

min _! > k, 2 > max 

v 2 |>0 v 2j U|i >0 u,, 



mm „; > k !( > max 

v, l >0 v,, u 21 >0 u 21 

All spectra uncertain. 

1 spectrum of ea comp. given 

by ext. k, other by int. k. 



v?|>0 v 2) 



Un>0 u,. 



mm _! > k„ > max 

v,j>0 v,, Uj,>0 Uj, 

All spectra uncertain, 

Botn spectra of ea component 

given by Intermediate k's. 
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Rank Estimation by Frequency Analysis 
of the Eigenvectors [9] 

The need to correctly determine the rank of a matrix for 
qualitative analysis has been demonstrated previously. Most 
rank estimation methods are statistical in nature and depend 
on the correct evaluation of the variance in the data. The 
method of rank analysis presented here uses Fourier Trans- 
form Image Analysis and is based on the assumption that the 
primary eigenvectors of an EEM (i.e. , those associated with 
the fluorescence) will contain information that is weighted 
toward the low frequency Fourier coefficients since spectra 
tend to be broad-banded. Similarly, the transforms of the 
secondary eigenvectors will have larger high frequency co- 
efficients. However, this manuscript is not a suitable 
medium for a comprehensive presentation of the theory of 
Fourier transforms. The reader is referred to any introduc- 
tory text on the subject of Fourier analysis [10]. 

The Fourier transform is fairly simple to implement since 
fast Fourier transform algorithms are routinely available for 
use even on small computers. Since any continuous function 
can be reproduced by addition of a series of sine and cosine 
functions with various frequencies, amplitudes, and phases, 
the forward Fourier transform is a method of determining 
these frequencies, amplitudes, and phases from the func- 
tion. The inverse transform reconstructs the time domain 
function from the frequencies, amplitudes, and phases. The 
discrete Fourier transform equation used to transform the 
eigenvector to the frequency domain is 

jV-1 
V(«) = UN 2 v(x)exj>-2TTi(xu)/N (lb) 

x=0 

where u is the frequency domain coordinate, V(u) is the 
Fourier transform of v(jc), and JV is the number of points in 
the discrete approximation of the function. The complex 
frequency domain function, V(«), is frequently represented 
by the Fourier spectrum which is given by 



\V(U)\ = [v(H) 2 real +V(H) 2 im J 1/2 



(2b) 



The area of the Fourier spectrum, 7\ is the sum of all the 
frequency coefficients and is given by 



N-l 

T = 2 i y (»)i 



(3b) 



K=0 



The segment of the Fourier spectrum bounded by +/— ulim 
is denoted A ulim which is defined as 



1 ulim 



w=wlim 

= 2 i y (")i 



(4b) 



«=Mlim 



which is simply the sum of the frequency coefficients in the 
section of the Fourier spectrum bound by ulim and -ulim. 
The relative importance of this frequency region in repro- 
ducing the time domain eigenvector, v(x), can be expressed 
by calculating the percent of T which lies in this frequency 
range. The parameter which represents the importance of 
the range from ulim to -ulim is called %A u]im and is calcu- 
lated from the expression 



%A uiim = A ufim /7"xl00 



(5b) 



When the value of ulim is well chosen, there is a marked 
drop in the %A ulim for secondary eigenvectors. Table 1 
shows a table comparing the rank estimation by frequency 
analysis to four statistical methods for matrices with ranks 
greater than or equal to 3. The frequency analysis method 
was the most accurate on the data tested here. 



Eigenvector Ratioing [11] 

This method was developed to resolve the spectra of 
compounds from the ellipticity matrices of two component 
mixtures. The presence of legitimate negative values in this 
matrix produced a need for a different algorithm to analyze 
this matrix despite its similarity to the EEM. The ellipticity 
matrix, F, is also bilinear and can be expressed as 



F=ST 



(lc) 



where the columns of the matrix S are the circular dichroism 
(CD) spectra, s k of the fluorophores in the sample, and the 
rows of the matrix T, symbolized by the vectors t k , are the 
emission spectra of those chiral fluorophores. 

The CD and emission eigenvectors of F are also an or- 
thonormal basis which span the matrix, F. The CD eigen- 
vectors, q t , may be represented as the columns of the matrix 
Q and the emission eigenvectors, p t , as the rows of the 
matrix P. Thus, F is also given by the equation 



F=OP 



(2c) 



assuming that the eigenvalues have been absorbed into one 
of the matrices. 

These eigenvectors must also be transformed to possible 
spectral vectors. However, the CD eigenvectors should not 
be transformed to a non-negative vector because the CD 
spectral vectors are not always non-negative. However, 
since the matrix F contains a finite amount of information, 
if the correct emission spectral vectors are found, the corre- 
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Table 1. EsliniLiletl numbers ol' components in text mixtures. 





No. 


of Components 






Kstimuled No. 


of Components 




Mixture No. 




K iiown 


REFAb 


Statl 


Stat2 




Si al .1 


Srat4 


1 




.i 


3* 


3* 


2 




3* 


4 


2 




3 


3* 


2 


2 




3* 


4 


3 




3 


3* 


2 


2 




2 


2 


4 




2 


2* 


1 


1 




2* 


2* 


5 




2 


2* 


1 


i 




2* 


3 


f> 




2 


T* 


2* 


2* 




?* 


.1 


7 




3 


-V 


3* 


i 




-p 


5 


H 




4 


4* 


3 


3 




3 


5 


<■} 




? 


5* 


3 


3 




3 


5 


ID 







5 


3 


3 




4 


(y' 



*" JeiMies a corrcrlly cslimnieil ruiik 
Slid I -- IJ^cJiviilue [X'lturhiilitm 
Sl;jl2 - Avcrii^L 1 emir nvlliiul 
Sul3 iv Chi •n]uaivJ iirsr 
SI.H4 - .%inn<h:U ermr lo-sc 



This is valid since 



F=Q'P' 



(9c) 



sponding CD spectral vectors are fixed and easily obtained. 
In this algorithm, the elements of the matrix, K, which 
transforms the emission eigenvectors to non-negative vec- 
tors are sought, Then, the possible CD spectral vectors are 
generated using the transformed emission vectors. The 
transformed emission vectors are the rows of the matrix P' 
which is given bv 



and 



P' = KP>0 



(3c) 



Far a two component mixture the values of the elements of 
K can be determined by expressing eq (3c) in terms of the 
elements of the matrices and solving the resulting expres- 
sions for the elements of K algebraically: 



P2/'=*22ft/ + *2lPli 2:0 



(4c) 
(5c) 



Again, assuming that k n -k 2 2 = \- tne elements of K are 

within the ranges defined by 



itiin 


Ph 5=A' r i& max 


~Pm 


(6c) 


Pv«> 


Pli Pli >0 


Pn 




min 


~Pn ^A-1,55 mas 


-Pn 


(7c) 


Pu <0 


Ph P\< >0 


Pu 





These ranges are generated because the sense of an inequal- 
ity changes when both sides of the inequality are divided by 
a negative number, 

The possible CD spectral vectors are generated by solving 
the following equation: 



Q=FP''(P'P'V 



(8c) 



I=P'P' r (P'P r ) 



(10c) 



This means that eq (8c) is another representation of the 
equation 



Q=Q'I 



(lie) 



Due to the weaker constraints on the transformation ele- 
ments in this algorithm, non-converging ranges are gener- 
ated for the values of the transformation elements. How- 
ever, it was found that if either of the components is a sole 
emitter in the monitored emission range, the value of the 
transformation element needed to transform an eigenvector 
to the spectra of the other component is usually given at an 
extreme of one of the ranges in eqs (6c) and (7c). This is 
because the values from the regions of sole emission best 
meet the criteria expressed in eqs (6c) and (7c). Therefore, 
except when the spectra are totally coincident at the base- 
lines at least one spectrum of each component is often re- 
trievable using this technique. Figure 2 shows an ambiguity 
table that was generated to illustrate the usefulness of this 
algorithm. 

This algorithm has also been tested on EEMs and it was 
found that this technique can be used with either set of 
eigenvectors: the emission eigenvectors (as it is with the 
ellipticity matrix) or the excitation eigenvectors. The results 
can be summarized hy an ambiguity table similar to the one 
in figure 2 which had excitation overlap on the vertical axis. 
The ambiguity table using only one axis is more ambiguous 
than the earlier table (fig. 1) using both axes. These results 
do not conflict with those found using eigenvector analysis 
on the EEM. "Multiplying" the two ambiguity tables gener- 
ated by this algorithm yields the table generated for eigen- 
vector analysis, verifying the validity of both methods. 
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Emission Overlap 



- Pz, " P2, 

mm — i. > k^ > max 

Pi,<0 P,, p,,>0 p,, 


- P?, - P?, 
mm 1L > k,, > max 

Pi,<0 P,, " P,,>0 p, 


- P?, ~ P?i 

mm — ii > k 12 > rfiax — — 
Pi,<0 p,, p,,>0 p„ 


- Pr " P21 
mln — X* > k, 2 > max 

P.,<0 Pi, Pi,>0 Pi, 


"Pi, -Pi, 

min ^ k 21 > max 


" p 'i , ~. " p " 
mm , , > k ?1 > max 


-Pi, , "Pi, 

mm > k S i > max 


" p, i -. , -^ "Pi' 
min ^ K 2i = max 


P;,<0 Pi, p 2 ,>0 Pj, 


P!,<0 P:, P;,>0 Pj, 


p 3 ,<0 p r , P?>0 P;., 


Pj,<0 p ; , p 2l >0 p 2 , 


All spectra uncertain 

Both spectra ot ea component 

given by extreme k's 


All spectra uncertain 

1 spectrum of ea. comp. given 

by ext. k. other by int. k. 


All spectra uncertain. 

1 spectrum of ea. comp. given 

by ext, k. other by int. k 


All spectra uncertain. 

Both spectra of ea. component 

given by intermediate k's. 



Figure 2— Ambiguity Table of Two-Component Ellipticity Matrix. 



Conclusions 



data using a diode array detector are bilinear in form. 
Clearly, there is a need for further development in this area. 



This manuscript has provided an overview of qualitative 
analysis techniques developed for matrix formatted fluores- 
cence data. Qualitative analysis of matrices was shown to 
generally consist of three basic procedures: rank estimation, 
determination of unknown component spectra, and screen- 
ing of expected compounds. The techniques outlined here 
only addressed the first two phases of the problem but they 
represent only a portion of the methods that have been 
developed to fill these requirements. The methods presented 
here successfully attack the stated problems within the 
framework of the limitations described. 

Rank estimation by frequency analysis can sometimes be 
more accurate than statistical methods for evaluation of 
spectral data. It would be useful to develop a criterion where 
the algorithm will automatically select a useful range for 
differentiating the secondary eigenvectors from the primary. 

In the present approach, eigenvector analysis and eigen- 
vector ratioing are limited to binary mixtures. Few real 
samples are binary. These methods must be extended to 
higher order mixtures. 

These techniques have been developed for use with fluo- 
rescence data, but are generally applicable to other forms of 
matrix formatted data. Some of the algorithms require that 
the data matrix be bilinear in form; however this is a charac- 
teristic of many types of data. For example, diode array 
detection of liquid chromatography and absorption kinetic 
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